Crossover Scaling in Dendritic Evolution at Low Undercooling 
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We examine scaling in two-dimensional simulations of dendritic growth at low undercooling, as 
well as in three-dimensional pivalic acid dendrites grown on NASA's USMP-4 Isothermal Dendritic 
Growth Experiment. We report new results on self-afhne evolution in both the experiments and 
simulations. We find that the time dependent scaling of our low undercooling simulations displays a 
cross-over scaling from a regime different than that characterizing Laplacian growth to steady-state 
growth. 
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Recent computational and experimental advances in 
dendritic growth offer a realistic prospect for a first 
principles understanding of solidification microstructure 
formation. Early experiments [Q,p| by Glicksman and 
coworkers on succinonitrile (SCN) provided the first 
benchmarks for theoretical models of dendritic growth. 
Comparison of experiments with theory has been diffi- 
cult, however, since experiments were influenced by con- 
vection effects and performed at low undercooling using 
materials with low anisotropy, parameters for which com- 
putation is difficult. Such calculations can nevertheless 
be performed in two-dimensions (2D) with state-of-the- 
art numerical methods combining so-called phase-field 
models |3l-|l3| and adaptive-mesh refinement |14-U(|. In 



the most recent round of experiments [^7 18 Glicksman 
and coworkers have reported observations on Pivalic Acid 
(PVA), whose higher anisotropy brings the benchmarks 
closer to the parameter range of theoretical computa- 
tions. 

Predicting dendritic growth theoretically has focused 
on the tip speed and shape in the steady state, simula- 
tions in 2D by Karma jl9| and subsequently ourselves |l5| ] 
have convincingly shown that the dynamically selected 
steady state is indeed the fastest of the discrete set of 
allowed needle crystals, as predicted by solvability the- 
ory pO|-p^|. However, at low undercoolings the diffusion 
length is so large that the time needed for each dendritic 
arm of a growing crystal to be in isolation from the others 
becomes much longer than any realistic simulation time 
[pTsfl . This regime, where dendrite arms cannot simply 
translate at a uniform speed because of their mutual in- 
teractions, was first systematically analyzed by Almgren 
et al. 1 25 1, who used solvability theory to explore the 



case when the temperature field strictly obeys Laplace's 



equation, with constant flux, at each time. They demon- 
strated that the dendrite tip position grows with the 3/5 
power of time, whereas the width grows with the 2/5 
power, results which were consistent with subsequent ex- 
periments in Hele-Shaw flow p6fl . 

In this Letter we explore dendritic growth dynamics 
at low undercooling, using the full diffusion equation 
dynamics. We find that the time-dependent evolution 
of 2D dendrite profiles is self-affine in time, generaliz- 
ing the results of Ref. j2?]] for the case of growth with 
a non-constant flux. Underlying this scaling behavior 
is a power law dependence on time of the dendrite tip 
position and maximum dendrite width. We find that 
scaling of these quantities displays a cross-over from a 
growth regime different from that of Hele-Shaw flow, to 
one characterized by steady-state tip growth. Meanwhile, 
comparison of our low undercooling simulations with mi- 
croscopic solvability theory gives good agreement for the 
value of the so-called stability parameter. We also exam- 
ine scaling in 3D dendrite data on pivalic acid obtained 
from NASA's USMP-4 Isothermal Dendritic Growth Ex- 
periment (IDGE), also finding sclf-affinc scaling in the 
global time-dependent PVA dendrite profiles. 

The simulated dendrites are modeled using the phase- 
field model employed in . Temperature T is rescaled 
to U = cp(T — Tm)/L, where cp is the specific heat 
at constant pressure, L is the latent heat of fusion and 
Tm is the melting temperature. The order parameter is 
defined by 4>, with cf> = 1 in the solid, and <f> = — 1 in the 
liquid, and the interface defined by 4> — 0. Id what follows 
time is rescaled by the time scale r c characterizing atomic 
movement in the interface, and length by the length scale 
W characterizing the width of the liquid-solid interface. 
The model is given by 
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where Z) = ar /W 2 , a is the thermal diffusivity, and 
A controls the coupling of U and 6. Anisotropy has 
been introduced in Eqs. (Q) by defining the width of 
the interface to be W(n) = W A(n) and the character- 
istic time by r(n) = T A 2 {n) @, where A(n) £ [0,1], 



and A(n) 
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The vector 
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n = {<P, x x + <t>,yy)/{<t>" x + tfy) 1 ^ 2 is the normal to the 
contours of <fi, and and <^ i2/ represent partial deriva- 
tives with respect to x and y. The constant e parame- 
terizes the deviation of W(n) from W . We expect the 
results to be similar for other definitions of anisotropy 
p8j. The parameters of Eqs. ([[]) are related to the ap- 
propriate Stefan problem using the relationships given in 
p9[ . In particular, W, r, A and D may be chosen to sim- 
ulate an arbitrary, anisotropic capillary length d(n) , and 
interface attachment coefficient fi{n), which we chosen 
here as j3 — 0, a limit appropriate for SCN and PVA. 

Simulated dendrites were computed by solving Eqs. (|l|) 
using the adaptive-grid method of Ref. |T^Jlq ]. Simu- 
lated dendrites were grown in a 2D quarter-infinite space 
using zero-flux boundary conditions along the sides of 
the system. Growth was initiated by a small quarter 
disk of radius R centered at the origin. The preferred 
growth directions are along the x and y axes, making 
these the directions of growth of dendrite branches. The 
order parameter is initially set to its equilibrium value 
4>o{x) = — tanh((|x| — Ji )/v2) along the interface. The 
initial temperature decays exponentially from U = 
at the interface to its far-field, undercooled value —A 
as f — > oo. Simulation data presented in this paper 
were obtained for three undercoolings: A = 0.25 and 
0.1 and 0.05. Details of these data are presented in Ta- 
ble 1. The two data sets for A = 0.1 correspond to 
different minimum grid spacings Aa; m ; n [|l5|,|l6|. Seed 
radii used in our simulations were R = 8.5, 15, 30, 30 for 
A = 0.25, 0.1(A), 0.1(B), 0.05, respectively. In all cases 
R is smaller than the thermal diffusion length by a factor 
of 20 or greater. 

The results of our low undercooling simulations are 
also contrasted here with new experimental data ob- 
tained from PVA dendrites. These experiments were 
performed by four of the authors (LaCombe, Lupulescu, 
Koss and Glicksman) during NASA's USMP-4 Isother- 
mal Dendritic Growth Experiment (IDGE). This exper- 
iment is described in detail elsewhere p8[. The IDGE 
experiment is designed to study dendrites grown under 
microgravity conditions, where transport in this partic- 
ular process is considered to be conduction-limited. The 
crystals are grown in an undercooled melt, controlled to 
within 0.001K. Growth is monitored thermometrically, 



while images are obtained from two perpendicular direc- 
tions using video and still cameras (electronic and film). 
Experimental results presented here were compared with 
four independent data subsets for dendrites grown at un- 
dercoolings of 0.58K, 0.63K, and 0.47K. We present the 
results from experiments corresponding to A = 0.052. 
These data were captured at times t\ — 42.48, t-i — 62.73, 
and £3 = 82.98 seconds after the dendrite was detected. 
The anisotropy for PVA was estimated at e pva = 0.025 
& 

We found the individual primary arms of our simu- 
lated dendrites to be sclf-affine, beyond some transient 
time, at all undercoolings examined. Figure 1 shows 
the (A-dependent) scaling profile for 2D dendrites grown 
at A = 0.05 and A = 0.25, respectively. The global 
scaling profile is obtained by scaling the x-direction by 
(x — Xb) I ^V max , where V max (t) is the distance from the 
tip Xtip(t) to the base Xb(t) of the dendrite arm, and the 
y-direction by J//V max (t), the maximum half- width of the 
lateral dimension of the primary dendrite arm. The tip 
and transverse directions were found to scale as V max ~ 
and V max - fT, where for A = 0.05, 0.1 A, 0AB and 
0.25, (3 = 0.0.73,0.73,0.78,0.97, and 7 = 0.43,0.43,0.45, 
and 0.55, respectively. For the A = 0.25 data, which 
at late times contained sidebranching induced by lattice- 
noise, we define V max (i) using the mean interface posi- 
tion, obtained by smoothing the data. This definition of 
the sidebranch envelope gives different results than using 
the maximum of the sidebranch envelope pO| , ^l[ . 

At low undercooling, long-lived transient interactions 
between neighboring primary dendrite arms causes their 
velocity and tip radius to deviate (within simulation time 
scales) from their steady-state values predicted by solv- 
ability theory ]T^| . However, we do find that the stability 
parameter a* — 2d D/VR 2 , where V and R are the time- 
dependent velocity and tip radius, agrees well with the 
value predicted by solvability theory. Figure 2 shows a* 
vs. time from our simulations at A = 0.25, 0.1 and 0.05. 
Error bars were estimated using AV, the fluctuations in 
velocity, and AR, deviations in radius of curvature. The 
radius was obtained by fitting to a second order polyno- 
mial near the tip. Deviations in the fit gave an estimate 
for AR. Data for A — 0.1 set B, omitted for clarity, con- 
verge to approximately the same a* as the A = 0.1 set 
A data, but display somewhat larger fluctuations around 
the mean, due to the larger grid spacing used. 

The time-dependent behavior of the tip position and 
lateral growth rate of our 2D dendrites are characterized 
by the scaling of V max (t) and Y max (t). Figure 3 shows 
V max and V max scaled onto respective crossover functions 
of the form 
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The parameters Lr> and td are effective diffusion length 
and time scales characterizing the intermediate regime, 
and are fit to give collapse of the A max and Y max data. 
The data for Fx(z) show a crossover scaling from fit 



to approximately Fx(z) 



-0.25 



at early times to 



Fx(z) 



in the steady-state regime. The cross- 



over in Fy(z) is given by Fy(z) 



-0.07 



at small z to 



Fy{z) 



v 0.05 



at large arguments of Fy(z). Exponent er- 



rors were approximately ±0.02, except for the A = 0.25 
data at late time, where they were ±0.05. These asymp- 
totic limits are demonstrated by the leveling off of Fx (x) 
and Fy(x) as x = t/ T D becomes large. 

Our simulations are in a regime where the dendrite 
is much smaller than the diffusion length. However, 
we do not observe the early-time scaling described by 
Almgren, et al [^5| since their calculations assume that 
dendrites are grown with a constant flux, whereas in 
our simulations the far field is diffusive with a spec- 
ified small undercooling. To illustrate this difference, 
let us assume that the rate of change of solid-fraction 
evolves as F ~ , whereby the solidified area A <~ t 1+ *. 
Since X max ~ f 3 , Y max ~ V , A ~ X max Y max , and so 
1 + / = j3 + 7. Since VR 2 = const, we obtain the scal- 
ing relation S(f,/3) = 4/ — 5/3 + 3 = 0. Our early-time 
exponents give 5(0.18,0.75) = -0.03 ± 0.1. In the late- 
time regime, 5(0.52,0.97) = 0.23 ± 0.22. The late-time 
error in 5 arises when we estimate 7 by smoothing the 
dendrite profiles for A = 0.25 (e = 0.05). For compari- 
son, we produced data for A = 0.25, e = 0.025 (not in 
Table 1), and which were free of spurious sidebranches, 
obtaining Y max ~ t 5 . 

Self-affine time-dependent scaling was also found in the 
mean dendrite profiles of the new 3D IDGE PVA data. 
Figure 4 shows the scaled PVA data for t — ^1,^2,^3- For 
comparison these data are superimposed on our 2D sim- 
ulation data for A = 0.05, e = 0.025. There is a slight 
asymmetry in the PVA data, likely due to interactions 
with other dendrite arms. For this reason, we scaled 
with respect to the top side the dendrite profile. Similar 
scaling was observed in all four IDGE data sets. The ex- 
perimental and simulated profiles show clear differences 
near the tip, as one would expect. Curiously, however, 
the profile shapes are in good agreement away from the 
tip. Similar results were also found in our A = 0.1 data. 
The reason for this is illustrated in Fig. 1, which shows 
that the difference in 2D dendrite profiles is small away 
from the tip. 

As plausible explanation for the apparent agreement 
between our low undercooling 2D simulations and our 
3D experimental data, we note that away from the tip, 
the diffusion field is more cylindrically symmetric than at 
the tip because the local diffusion length is larger. Thus 
diffusion of heat away from the interface is better approx- 
imated by the 2D diffusion equation. We hope to examine 
this idea critically, as well as to accurately determine the 
experimental scaling behavior in future publications. 
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TABLE I. Parameters for simulated dendrites. The time 
t* = 255622.4 
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FIG. 3. Crossover scaling functions describing lateral width 
of simulated dendrite arm Y max and tip-to-base distance 
A max , for A = 0.25, 0.1 (sets A and B) and 0.05. 
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FIG. 1. Comparison of scaled dendrite profiles for A = 0.05 
and A = 0.25. For A = 0.25, nine times are plotted, spaced 
between 28643 < t < 66083. For A = 0.05, there are six times 
in the range 222022 < t < 279622. 
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FIG. 4. Scaled USMP-4 PVA dendrites grown at 
A = 0.052. 2D simulation data is superimposed for A = 0.05, 
e = 0.025. Comparison of 2D and 3D data made merely to 
illustrate self-affinity in both 2D and 3D dendrites. 
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FIG. 2. Simulation data of a* vs. time for A = 0.25,0.1 
(set A) and 0.05. For clarity, the A = 0.1 and 0.25 data have 
been shifted along the y-axis by 0.04 and 0.02, respectively. 
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